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Superconducting Order Parameter in Bi-Layer Cuprates: Occurrence of n Phase 

Shifts in Corner Junctions 
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We study the order parameter symmetry in bi-layer cuprates such as YBaCuO, where interesting 
7r phase shifts have been observed in Josephson junctions. Taking models which represent the mea- 
sured spin fluctuation spectra of this cuprate, as well as more general models of Coulomb correlation 
effects, we classify the allowed symmetries and determine their associated physical properties, ir 
phase shifts are shown to be a general consequence of repulsive interactions, independent of whether 
a magnetic mechanism is operative. While it is known to occur in d-states, this behavior can also 
be associated with (orthorhombic) s-symmetry when the two sub-band gaps have opposite phase. 
Implications for the magnitude of T c are discussed. 

PACS numbers: 74.72. -h , 74.20.Mn, 74.50. +r, 74.62. -c 



The observation in YBCO of unusual Josephson junc- 
tion behavioitTu is one of the most important experimen- 
tal results to emerge from the cuprate literature in recent 
years. Here in a Josephson SQUID experiment the two 
junctions are configured so that their normals lie along 
the a and b axes of the Cu02 plane. This measurement 
has been widely interpreted as support for a d-symmetry 
of the order parameter, as well as for a magnetic mech- 
anism for superconductivity. In this paper we show that 
both of these inferences may be inappropriate. For nota- 
tional precision, throughout this paper we use the terms 
s- ( or d- ) symmetry to correspond to states which have 
the same (or opposite) sign under a 7r/2 rotation. Thus 
the d-states under consideration are not necessarily of 
the specific d x 2_ y 2 form. 

The gap equation for bi-layer systems has been stud- 
ied earlier in the xantext of a magnetic mechanism for 
superconductivity™! There it was observed that the 
d-symmetric state of the single layer problem, is trans- 
formed to a pair of in-phase d states on each of the two 
sub-bands, and that these compete with a pair of out- 
of-phase s-states. Here we take the problem to a greater 
level of generality, establishing that this situation persists 
for a wide class of repulsive interactions, which are un- 
related to the antiferromagnetic spin fluctuation picture. 
Alternate classes of the order parameter symmetry are 
also generated. These correspond to in-phase s-states 
and out-of-phase d-states. We establish how the rela- 
tive stability of the two competing states is affected by 
changes in bandstructure, orthorhombicity, and hole fill- 
ing. 

It should not be surprising that d-states have a more 
general origin beyond the antiferromagnetic spin ex- 
change models. In a one layer cuprate, the lattice symme- 
try requires that all gap states are either even (s-) or odd 
(d-) under a 7r/2 rotation. In bi-layer systems, these one 
layer states generalize naturally to a pair of even or odd, 
in-phase or out-of-phase states, associated with each of 
the two sub-bands. Thus, as one of two allowed states, 
d-symmetry should be widespread, and independent of 



the microscopic details of the model. 

In the presence of both intra- and inter-layer interac- 
tions (Vj| and Vj_), the weak coupling BCS gap equation 
becomes a set of coupled equations for the gaps on each of 
the sub-bands. It is simpler to write the gap equations in 
terms of the two sub-band gaps A + , A_ rather than the 
intra-layer (Ay) and inter-layer (A^) components. These 
are related via the unitary transformation which diago- 
nalizes the Hamiltonian. In this model t± is the matrix 
element for hopping between layers; t and t' refer to the 
first and second nearest neighbor in-plane hopping which 
may contain orthorhombic effects. On site Coulomb ef- 
fects U are assumed to enter via a renormalization of the 
bandstructure parameters as shown by Si et aH. 

Following the usual procedure!!!, the gap equations be- 
come 
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where the superconducting quasi-particle energies are 
E± = 



. + A j_ , where 

e± = — 2t[cos(q x a) + cos(q y b)] 

+ 4t' cos(g x o) cos(q y b) ± t_ 



(2) 



It follows that each of the two sub-band gaps can be 
written in terms of the parallel and perpendicular com- 
ponents as 



A,, = (A+ + A_)/2, A x = (A + - A_)/2 



(3) 
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In the case of a magnetic pairing mechanism, the two 
interactions are related to components of the dynamical 
spin susceptibility. This susceptibilitit-has been calcu- 
lated for the bi-layer cuprate YBCOQa, in the strong 
U limit. For realistic Fermi surface shapes, and mod- 
erate in-plane and out of plane exchange interactions, 
the results .axe in reasonable agreement with neutron 
experimental!]. Because the in-plane magnetism is not 
independent of inter-plane effects any proper treatment 
of spin fluctuation induced superconductivity should in- 
corporate both components. A reasonable approximation 
of both tho-theoretical resultsu and the experimental neu- 
tron dataE9 is to take 

Vj|=$jfx||(q-q') (4a) 

Vx =sixx(q-q / ) 

= -£Lxu(q-q') (4b) 

with 

X|1 {kx ' ky) = [l + J { C osk x a + C osk y b)Y (4C) 

and J « 0.3. Here we have absorbed overall coefficients 
into the coupling constant prefactors gn and gx- This 
model is similar to that used in Ref.5, except that we have 
assumed an arbitrary relation between the magnitudes of 
two superconducting coupling constants, which are taken 
to be the same in Ref.5. Moreover, we fit the dynamical 
susceptibility to neutron, rather than NMR dataEH. 

It may be IMted from Eq. ([k]) (as well as experimental 
neutron datalij) that the antiferromagnetic fluctuations 
show up as a weak peak around (7r/a,7r/a) with short 
coherence lengths, suggesting that the system is far from 
any real instability; thus high T c is difficult to explainEil. 
Consequently we explore more general mechanisms by 
extending Eqs. (|Ia| ) and (^) to the case where the over- 
all signs are unconstrained and the q — q' peaks occur at 
arbitrary wave- vector with arbitrary peak width J Q . We 
define 

y _ ^A, Ag . . 

l± [l-J (cos(k x ±Q x )+cos{k y ±Q y ))}2 {> 

where all signs are summed over. We divide our numer- 
ical analysis into four distinct cases in which either the 
inter-layer or the intra-layer interaction dominates and 
in which the respective interaction is repulsive (A > 0) 
or attractive (A < 0).The model is viewed as a gen- 
eral representation of pairing mechanisms of the elec- 
tronic and phononic variety. Since the former usually 
derives from a generalized susceptiblity, and the latter 
from a phonon propagator, there are no sign changes as 
a function of momentum transfer. This poses impor- 
tant constraints on the allowed superconducting states. 
Our resultslI3 are summarized in Figures 1 and 2 for the 
case of characteristic wave- vectors q— q' along (%/a, ir /a) 
and along (jr/a,Tr/2a). The latter wave-vector illustrates 



the behavior away from the antiferromagnetic instabil- 
ity model, in order to show the generality of our results. 
Plotted in these figures are the form of the gap functions 
in the two regimes. The figures on the left (right) in each 
box correspond to intra-layer (inter-layer) dominated be- 
havior. Our conclusions from both Figures 1 and 2 may 
be succinctly summarized. We find that d-symmetry is 
associated with repulsive and s-symmetry with attrac- 
tive interactions. In phase gap behavior occurs when the 
intra-layer interaction is the larger; out of phase behavior 
arises in the opposite case. This phase dependence can 
be deduced from Eqs.(|J). In the case of dominant inter- 
layer effects, |AjJ > |A|||. This will occur when A + and 
A_ have opposite signs. 

To establish the generality of these results, we have 
varied the Fermi surface shape ( via the ratio of t'/t ), the 
position of the Fermi energy or hole filling and the width 
J of the peak structure. These variations introduce only 
quantitative but not qualitative changes in the physical 
picture shown in the two figures. 

It is important to note from Figs.[l](a)-(c) and ||(a)-(c) 
that the out of phase s- and d-wave states wjJl|-exhibit 
7r phase shifts in a corner SQUID experimentErta. This 
corresponds to a change in sign of the " sum " order pa- 
rameter A|| of Eqs.(||) upon varying from to tt/2. While 
not a general feature of all solutions, its presence requires 
(a modest amount of ) orthorhombicity. The observation 
of 7r phase shifts in bi-layer cuprates is thus not as strong 
a constraint on the order parameter symmetry as in one 
layer materials. All orthorhombic states which exhibit 
these 7r phase shifts will also show finite c-axis tunnel- 
ing in untwinned crystals!!! However, twinning effects 
(if they average fully over the a and b-axes), will lead to 
a cancelation of Josephson coupling, whenever the corner 
SQUID experiment has the observed it phase shift. 

It should be stressed that the out of phase s-states have 
the additional advantage as a candidate gap symmetry, 
over d-states (in or out of phase), of being relatively in- 
sensitive to impurity effects. In addition^— this state can 
be compatible with neutron experimentstLO, which show 
no nodal signature. In contrast to experiment, because 
nodes are not present, at least in the. clean limit, power 
laws in thermodynamical propertiealj are not expected. 

We have searched for nodal behavior in these s-wave 
states with spaie care, since there is recent photoemis- 
sion evidence&3 to suggest that they may exist in bi-layer 
BISCO. Several observations are important to note in 
this context. (1) As the inter-layer hopping t± becomes 
small the magnitudes of the two gaps become equal and 
they are less able to respond to orthorhomicity by pro- 
ducing gap anisotropy. Consequently BISCO 2212, which 
is believed to have a very small tx, would be unlikely to 
exhibit nodal s-wave behavior. (2) The eight node s-stal£ 
which has been conjectured as a candidate for BISCOE£l 
appears quite generally as a meta-stable state whose so- 
lution has a lower T c than the nodeless s-wave or (four 
node) d-wave symmetric gap. Within the manifold of 
meta-stable states, the more nodes, the 
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FIG. 1. Superconducting gap for interactions peaked at q — q' = (7r,7r), for the case of attractive and repulsive intra-layer 
( Aa) and inter-layer (As) interactions. Figures on left (right) are for intra- (inter-) layer dominated regimes. 



FIG. 2. Superconducting gap for interactions peaked at q — q' = (7r,7r/2). 
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lower is the T c . (3) We have studied solutions to a 
separable pairing potential model in which the suscep- 
tibility of Eq.(f4c|) is replaced by a product of cosine 
terms: cosLQ^a) cos(q' x a) + cos(q y b) cos(q' y b) , since it was 
speculated!!^ that this potential would give rise to an 8 
node s-state. We find that d-wave states arise naturally 
in this model as well, and they are always more stable 
than s-wave states. 



FIG. 3. Effect of inter-layer correlation on the super- 
conducting transition temperature T c . (A) indicate the 
(d,d) states and (O) represent the (s,—s) states. We con- 
sider different shapes of Fermi surface and Fermi energies 
- orthorhombic lattice with next-nearest-neighbor hopping, 
Ef = —0.2 (dotted line), Ef = —0.1 (solid line); and tetrag- 
onal lattice with t' = 0, Ef ~ —0.6 (dot-dashed line), 
E F = -0.4 (dashed line). 

In Fig.|| is plotted the dependence of the transition 
temperature on the inter-plane coupling constant T c for 
the magnetic pairing r«iodel of Figj^(a), for two different 
Fermi surface shapesQ corresponding to YBaCuO (with 
two hole concentrations as well ) and LaSrCuO (t' = 
). Here, because we use a weak coupling approxima- 
tion, the absolute values of T c are not meaningful. How- 
ever, the relative changes with different parameteriza- 
tions are expected to be accurately captured. The cir- 
cles correspond to the out-of-phase s-states and triangles 
to in-phase d-states. As expected, for sufficiently small 
inter-plane coupling the d-wave state is the more stable; 
however depending on the bandstructure and Fermi sur- 
face shape, this state may be readily de-stabilized to a 
pair of s states by the introduction of a very small inter- 
plane interaction. This reflects the general result that 
(nodeless) s-wave states are able to take better advan- 
tage of the superconducting interaction than can d-wave 
states, which require a cancellation of positive and nega- 
tive terms to satisfy the gap equation. Note that in the 
LaSrCuO model the T c 's are generally higher as a conse- 
quence the better Fermi surface nesting along the direc- 
tion of the wave-vector (tt, 7r). For each parameter set, 
the various curves tend to coalesce at higher As, where 
the (s-like) states are found to be more isotropic. This is 
a consequence of the fact that these isotropic states are 
not able to utilize the Van Hove singularity effect, which 
is relatively more important for d x 2_ y 2 state. 



In summary, by solving the gap equation for bi-layer 
models with general repulsive interactions, we find that 
& x 2_ y 2 states arise quite generally and are not uniquely 
associated with wave- vector structure along the antiferro- 
magnetic direction. Moreover, we have established that 7r 
phase shift behavior, which is often cited as the strongest 
evidence for d-wave pairing can also be associated with 
(orthorhombic) s-symmetry when the two sub-band gaps 
have opposite phase. This state has some advantages 
over d-states in large part because of the relatively small 
sensitivity of T c to non-magnetic impurities. An impor- 
tant conclusion from our analysis is that there are al- 
ways competing states in bi-layer systems, and that the 
order parameter symmetry would be expected to vary 
from cuprate to cuprate as well as within a given cuprate 
class at different hole concentrations. One can conclude 
that the Josephson junction data, in particular, provide 
strong evidence for superconductivity mediated by some 
form of repulsive interaction. On the other hand, these 
collected observations (in bi-layer cuprates) weaken the 
often cited support for theories of spin fluctuation medi- 
ated superconductivity. 
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Note Added: 

After this manuscript was submitted we received a 
preprint from K. Kuboki and P. A. Lee, in which an 
RVB description of bi-layer superconductivity was used 
to infer a spontaneous breaking of tetragonal symmetry. 
This s, d mixing will not occur in the present model, as a 
result of the free energy form which contains a quadratic, 
rather than quartic, mixing of the in- and out- of- plane 
gaps. 
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